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Abstract 

Given two absolutely continuous nonnegative independent random variables, we de¬ 
fine the reversed relevation transform as dual to the relevation transform. We first apply 
such transforms to the lifetimes of the components of parallel and series systems un¬ 
der suitably proportionality assumptions on the hazards rates. Furthermore, we prove 
that the (reversed) relevation transform is commutative if and only if the proportional 
(reversed) hazard rate model holds. By repeated application of the reversed relevation 
transform we construct a decreasing sequence of random variables which leads to new 
weighted probability densities. We obtain various relations involving ageing notions and 
stochastic orders. We also exploit the connection of such a sequence to the cumulative 
entropy and to an operator that is dual to the Dickson-Hipp operator. Iterative formu¬ 
lae for computing the mean and the cumulative entropy of the random variables of the 
sequence are finally investigated. 

Keywords: Relevation transform; reversed relevation transform; proportional hazards 
rate model; proportional reversed hazards rate model; weighted cumulative distribution; 
cumulative entropy. 
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1 Introduction 

In this paper we introduce the reversed relevation transform and study some properties of a 
new weighted cumulative distribution function and its connection with the cumulative en¬ 
tropy. The considered stochastic model is dual to the weighted tail distribution studied by 
Kapodistria and Psarrakos |18| . Specifically, we construct a sequence of stochastically de¬ 
creasing random variables {X n ,n > 1}. In this sequence, X\ is nonnegative and absolutely 
continuous, and the (n+l)th random variable of the sequence is inductively defined through 
the following relation: [X n+ \ \ X n = t] = [ X n \ X n < t], n = 1, 2, ■ ■ ■ , for t > 0. Here, as 
usual, [X | B] denotes a random variable having the same distribution of X conditional on 
B , and = denotes equality in distribution. 

Roughly speaking, {X n ,n > 1} is suitable to describe an iterative process involving a 
sequence of tasks, where X n is the random time required to perform the nth task. For 
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instance, we refer to a training procedure based on iterative learning or a working system 
based on replacements or repairs of failed items. Starting from the sequence {X n ,n > 1}, we 
derive some properties of a new weighted cumulative distribution via stochastic orders, and 
its connection with covariance and cumulative entropy. Our investigation is also devoted to 
disclosing iterative rules that allow us to compute the mean and the cumulative entropy of 
the random variables of the sequence, whose computational efficiency is illustrated by some 
numerical examples. 

We recall that [18] constructed a stochastically increasing sequence of random variables, 
whose iterative rule involved the residual lifetime of X n , i.e. 

[X n+ i - X n I X n = t] = [X n - 1 1 X n > t\, n > 1, t > 0. 

They obtained some results on this sequence and its connections to the cumulative residual 
entropy. Their process may describe the successive failures of a component, which, on 
failure, is replaced by a component of equal age, but the lifetime distribution of the nth 
component is assumed to be identical to the distribution of the time until the nth failure; see 
[18] for more details. A simpler case was studied by Baxter 0, who considered a stochastic 
process generated by the successive failures of a component which on failure is replaced by 
a component of equal age. 

The paper is organized as follows. In Section 2 we present the reversed relevation 
transform with a preliminary result based on the usual stochastic order and an application 
to parallel systems involving the proportional reversed hazards rate model (PRHRM). Some 
dual results for the relevation transform are then provided. We also address the problem of 
determining necessary and sufficient conditions such that the reversed relevation transform 
and the relevation transform are commutative. The new weighted distributions and their 
characteristics based on stochastic orders and ageing properties are discussed in Section 3, 
where the new notion of decreasing reversed hazard rate (DRHR) in a length-biased sense is 
also considered. In Section 4 various kinds of entropy such as Shannon entropy, cumulative 
entropy and dynamic cumulative entropy are examined. Then the connections between 
the earlier mentioned entropies and several functions of the given sequence are discussed. 
Specifically, we also obtain some iterative results for the involved quantities. These include a 
new probabilistic meaning for the cumulative entropy, which can be expressed as a difference 
of means of consecutive random variables of the considered sequence. Finally, in Section 5 
we define an integral operator and we discuss its properties related to the previous results. 
Also, various numerical examples are presented to shed further light on the characteristics 
of the studied sequence. 

2 Background and preliminary results 

Let X be an absolutely continuous nonnegative random variable with probability density 
function (PDF) f(t ), cumulative distribution function (CDF) F(t) = P(X < t ), and survival 
function F(t) = 1 — F(t), so that X may be viewed as the random lifetime of a system or 
a component or a living organism. Assume that F(t) > 0 for all t > 0. We recall that the 
reversed hazard rate function of X is defined by 

r(t) = -tlogm = M t> 0. 

There are several papers on applications of reversed hazard rate function in the literature, 
see e.g. Block et al. [4j, Di Crescenzo D33, Gupta and Nanda m, Kijima and Ohnishi [19l . 
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and the references therein. For the random lifetime X, we define Xm = [t — X\X < t], t > 0, 
which is termed the inactivity time. Indeed, Xw describes the length of the time interval 
occurring between the failure time X and an inspection time t, given that at time t the 
system has been found failed. For t > 0, the mean inactivity time of X is given by 


jl{t) = E[t -X\ X <t] 



( 1 ) 


It is known that the CDF of the past lifetime [X \ X < t], t > 0, is given by 


P(X < x\X < t) 


F(x) 

m 

i 


0 < x < t 

X > t, 


so that the PDF of the past lifetime is f{x)/F(t) for all 0 < x < t. 

Hereafter, we consider two nonnegative absolutely continuous and independent random 
variables X and Y with the CDFs F(t) and G(t), respectively. 


Definition 1 The reversed relevation transform of X and Y is defined by: 


GftF(x ) — |-^yl{0<a;<t} + l{x>t} jdG(t) 

= G(x) + F(x) J -^yd G(t), x > 0, (2) 

where 1,4 is the indicator function of the set A, i.e. 1 a(x) = 1 if x € A, and 1 a(x) = 0 if 
x € A c . 


Generally, the inactivity time of X at a random time Y. denoted by X[yj, is dehned as 

X[Y] = [Y — X \X < Y] (see, e.g. [18]). Moreover, let X\Y] = [X \ X < Y] denote the total 
time of X given that it is less than an independent random inspection time Y. Therefore, 
the CDF of X[Y] is given by 

P(X[T] < x) = G#F(x), (3) 

where the symbol ft is defined in (|2|). If random variables X and Y are independent and 
identically distributed (i.i.d.), then 

P(X[T] <x) = FftF(x) = F(x)[l + T(x)], x > 0, 


where 


r 

J X 


T(x) = — logF(x) = / r(u)du, x > 0, 


denotes the cumulative reversed hazard rate function of X\ see, e.g. [HO]. 


( 4 ) 


Example 1 (a) Let X and Y be independent nonnegative random variables having the cdfs 
F(x) = exp (—ax -7 ), x > 0 , and G(x) = exp (—cx -7 ), x > 0 , respectively, with a, c, 7 > 0 . 
From © and ([3]), we have 


F(X[Y] < x) = GftF{x) 


-(cexp(— ax 7 ) — aexp(—cx 7 )), 

c — a 

(l + cx -7 ) exp(—cx -7 ), 


x > 0 , a ft c, 
x > 0 , a = c. 
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(b) If F(x) = e a /( eX - 1 ), x > 0, and G(x) = e c ^ eX l \ x > 0, with a, c > 0, then 


x > 0, o^c, 
x > 0, a = c. 

Ageing notions and stochastic orders have many applications in various areas of science 
such as reliability and survival analysis, economics, insurance, actuarial and management 
sciences, and coding theory; see Shaked and Shanthikumar [32] for a detailed account. In 
the following, we review some notions that are used in the sequel. Note that here and 
throughout this paper, the terms ‘increasing’ and ‘decreasing’ are used in a nonstrict sense, 
and R denotes the set of real numbers. 

Definition 2 If A is an absolutely continuous random variable with support (l x ,u x ), CDF 
F, PDF / and reversed hazard rate function r(f) = f(t)/F(t ), then 

• X is said to have the increasing likelihood ratio (ILR) property if f(x) is log-concave 
or, equivalently, the function f'(x)/f(x) is decreasing in x £ {l X ,u x )', 

• X is said to have the decreasing likelihood ratio (DLR) property if f(x) is log-convex 
or, equivalently, the function f(x)/f(x) is increasing in x £ (l x ,u x ); 

• X has the DRHR if r(t) is decreasing in t £ (l x ,u x ) or, equivalently, T(x) = 
— log F (x) is convex. 

Moreover, if Y is an absolutely continuous random variable with support (l Y ,u Y ), CDF G 
and PDF g, then 

• X is smaller than Y in the usual stochastic order (denoted by X < st Y) if F(t) < G(t), 
for all t £ R, or equivalently F(t) > G(t), for all 1 £ K; 

• X is smaller than Y in the likelihood ratio order (denoted by X < lr Y) if f(x)g(y) > 
f{y)g{x) for all x < y, with x,y £ ( lx,u x ) U {l Yl u Y )] 

• X is smaller than Y in the up-shifted likelihood ratio order (denoted by X <i r ^ Y) 
if X — x <i r Y for all x > 0 or, equivalently, for each x > 0 we have g(t)/f(t + x) 
is increasing in t £ {l x — x,u x — x) U ( l Y ,u Y ), where a/0 is taken to be equal to oo 
whenever a > 0. 

Proposition 1 If X and Y are independent nonnegative random variables, then 

X\Y] < st min{A, Y }. 

Proof. Due to © and ©, for x > 0, we have 

P(X[Y]<x) = G(x) + F(x)j^^ r) dG(t) 

POO 

> G(x) + F(x) / dG(t) 

J X 

= G(x) + F(x) - F(x)G(x). 

The proof thus follows recalling that the last term is the CDF of min{X, Y}. 


P (X[Y\ <x) = G#F(x ) = < 


1 


c — a 
1 + 


( Ce -a/( e ^-l)_ a e-c/( e ^-l)), 


e x - 1 


^-c/(e x - 1) 
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Let us now consider a stochastic model that extends both cases treated in Example |T] 
Let X and Y be absolutely continuous nonnegative random variables with CDFs F(x) and 
G(x), and reversed hazard rate functions tx(x) and Ty(x), respectively. These variables 
satisfy the PRHRM with proportionality constant 9 > 0 if Ty(x) = 9tx(x) for all x > 0 or, 
equivalently, if 

G{x) = [F{x)] e , x > 0. (5) 

The parent distribution function can be expressed as F(x) = e~ T ^ x \ x > 0, where T(x) 
is defined in (f4j) . The model (JSJ) was first proposed by Lehman |23j in contrast to the 
proportional hazard rate model. It is more flexible to accommodate both monotonic as well 
as nonmonotonic failure rates even though the baseline failure rate is monotonic. For more 
details on the applications and properties of the proportional hazard rate model see e.g. 
m, na, m, m and m, among others. 

Proposition 2 Under the PRHRM (0), we have 

G#F(x) = (0e~ T W - e-^1) , x > 0 , (6) 

for 9 > 0, 9^1. 

Proof. Under the PRHRM ([2]), we can verify that 

rm iG{t) =^ to --**-™). 

This identity and d^J) thus yield □ 


Note that the assumptions of Proposition [2] are satisfied by the cases (a) and (b) shown 
in Example [1] for a ^ c. 


Example 2 Let X m:m = max{Xi,..., X m j be the lifetime of the parallel system consist¬ 
ing of m components with absolutely continuous i.i.d. lifetimes Xj., X m . having the 
common cdf F. Moreover, suppose that X t . i = 1,..., m, and Y satisfy the PRHRM with 
proportionality constant 6 , as in ([5]). Then after some calculations, from Q we can express 
the CDF of X m:m [Y] as the following generalized mixture: 


G#F m:m (x ) 


6F m {x) — mF e {x) 
9 — m 


x > 0, 


for 9 ^ rn. □ 

In the forthcoming theorem we investigate the commutative property for the reversed 
relevation transform. 

is commutative if and only if 

then Q yields 

x > 0, 


Theorem 1 The reversed relevation transform of X and Y 
X and Y satisfy the PRHRM. 

Proof. If X and Y satisfy the PRHRM as specified in (J5|). 

6F{x) - G{x) 


G#F(x) = F#G(x) = 


0-1 
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for 9 > 0, 9 7 ^ 1, and the desired result follows. To prove the converse, we assume that, for 
all x > 0 , 

/‘oo i roo i 

G(x) + F{x) J G{t) = Fix) + G{x) J ^yd Fit). (7) 

Differentiating both sides of (|7|) . we have 

/‘OO 1 f‘ OO | 

/w x m iG(t)= 9 (x) i* W) AF[t) ' (8) 

where / and g denote the PDFs of X and Y, respectively. Again, differentiating both sides 
of JSJ), we obtain 



fjx)gjx) 

Fix) 



fjx)gjx) 

G(x) 


From (ED-© and some algebraic simplification, we obtain 


(9) 


fix) _ fjx) = g’jx) _ gjx) 
fix) Fix) gix) Gixf 


Ji ln IM. = fL m 9jx) 

dx Fix) dx Gix )’ 

Integration on both sides yields 


x > 0 . 


ln 


fix) 

Fix) 


= ln 


gjx) 

Gix) 


+ constant, 


x > 0 , 


or 


fix) „ gjx) 
Fix) Gix) ’ 


x > 0 , 


where 9 is a positive constant. Thus, we obtain G(x) = [Fix)] 6 , x > 0, which completes 
the proof. 


Hereafter, we analyse some results that are dual to those given above. Let X and Y be 
absolutely continuous nonnegative random variables with survival functions Fix) and Gix), 
and hazard rate functions hxix) = — (d/dx) log Fix) and hy(x) = — (d/da:) log G(x), re¬ 
spectively. Then X and Y satisfy the proportional hazards rate model with proportionality 
constant 9 > 0, if hyix) = 9hxix) for all x > 0. This is equivalent to the model 

Gix) = [Fix)] 6 , 9 > 0, (10) 

where Fix) = e~^ x \ x > 0, is the parent survival function and A (a;) = — log-F(x), denotes 
the hazard function; see, e.g. [15]. Let X(Y) denote the total time of X given that it exceeds 
an independent random inspection time Y, i.e. X(Y) = [A|A > Y]. Then we have 

P(X(Y) > x) = G#F(x) 

= Gix) + Fix) J* J^dGit), x > 0, (11) 

where the symbol ff denotes the relevation transform introduced by Krakowski m- Equa¬ 
tion ED was discussed in [18]; see also [ 6 ] and m and the references therein. 

The following result is analogous to Proposition [T] and, thus, the proof is omitted. 
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Proposition 3 If X and Y are independent nonnegative random variables, then 

X(Y) > st m&x{X,Y}. 

Let us now see the analogous of Proposition [2J 
Proposition 4 Under the proportional hazards rate model m, we have 

G#F(x) = j£ (de~ A ( x) - e~ eA (x) ) , x > 0, 

ford >0,6^ 1 . 

Proof. Since 

)• x>0 ' 

the proof follows from dill) . □ 


Example 3 Suppose that Xi- m = min{Ai,..., X n } is the lifetime of the series system con¬ 
sisting of m components with absolutely continuous i.i.d. lifetimes X\,... , X rn , having the 
common CDF F. Also, suppose that A,., i = 1,... ,m, and Y satisfy the proportional haz¬ 
ard rate model with proportionality constant 0 , as in ( 1101 ) . By means of some calculations, 
from m we obtain the survival function of X\ :m (Y), which is expressed as a generalized 
mixture: 


G#F 1:m (x) 


8F m (x ) -mF e {x) 
8 — m 


x > 0 , 


for 8 yl m. 


□ 


The relevation transform is not always commutative. Indeed, in the following theorem 
we give a necessary and sufficient condition leading to such a property. Being similar to 
Theorem |T] we provide only a sketch of the proof. 

Theorem 2 The relevation transform of X and Y is commutative if and only if X and Y 
satisfy the proportional hazard rate model. 


Proof. Let X and Y satisfy the proportional hazard rate model as in (HOD . From (ED, 
we thus have 

GjfF(x) = FffG(x) = 6F ^~^ X \ x>0 ^ 

for 8 > 0, 8 7 ^ 1, and then the relevation transform is commutative. To prove the converse, 
we assume that for all x > 0 

G(x) + F(x) £ -^yd G(t) = F(x) + G(x) £ ^d F(t). 

Differentiating both sides and after some calculations, we obtain 


so that 


f'{x) f{x) _ g'(x) g[x) 
f(x)^F(x) g(x)^G{xY > ’ 


4-hiS = T ln »hl, x>0 . 


~— ill “ = —— ill -=- , 

dx F(x) dx G(x) 

Such a relation implies that G(x) = [F(x)] 61 , x > 0, for 8 > 0, thus completing the proof. 
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The characterization of distributions based on the relevation transform has been the 
object of various investigations; see, e.g. Lau and Prakasa Rao m and m ■ We point out 
that [14, Theorem 9] states that the relevation transform of two i.i.d. nonnegative continuous 
random variables is identically distributed to their convolution, i.e. F#F(x) = F *F{x) for 
all x > 0, if and only if they have exponential distribution. However, due to Proposition [1] 
a similar result cannot hold for the reversed relevation transform. 


3 Sequence of weighted distributions 


Let X be an absolutely continuous nonnegative random variable with PDF f(x) and CDF 
F(x). Based on X, we construct a sequence of random variables {X n ,n > 1} as 

X! = X, [X n+1 \X n = t\ = [X n \X n <t], n> 1, (12) 

or equivalently 

X!=X, X n+1 = [X n \X n <X’ n \, n> 1, (13) 

where X' n is an independent copy of X n . It is easy to show that the corresponding CDFs 
F n (x) = P (X n < x) are given as, for all x > 0, 

Fi(x) = F(x), F n+ i(x) = F n #F n (x ) = F n (x)[ 1 +T n (x)], n > 1, (14) 


where 


T n (x) = -logF n (x) = / T n (u)du, x > 0, 


(15) 


denotes the cumulative hazard function of X n , and 


T n {u) = u> 0 , (16) 

F n (u) 

is the reversed hazard rate of X n . From (ED, we can see that the corresponding densities 
are given by, for x > 0 , 

fl(x) = f(x), f n +l(x) = T n (x)f n (x), tl > 1. (17) 

Due to ED, for x > 0 , we have 

T n+1 (x) = — log F n+ i(x) = T n (x) - log(l + T n (x)), n> 1, 
and, thus, from (ED, we obtain 


n 


fn+l{x) = T n (x)f n (x ) = T n (x)T n _x(x)f n -l(x) = ■■■ = ]^[ Tj(x)/(x), n > 1. (18) 


i=1 


From ED, we see that fn+i(x) is a sequence of weighted PDFs. We recall that, given an 
absolutely continuous random variable X having density / and a nonnegative real function 
w, the associated weighted random variable X w has the PDF 


= w(x)f(x) 
- EHX)] 


x G R, 


provided that 0 < E[tu(X)] < oo. See [25], [29], and [ 2 ] for some recent papers on weighted 
distributions. 




The sequence of random variables {X n , n > 1} is suitable to describe an iterative process, 
where X n denotes the random time required to perform a task at the nth stage. For instance, 
consider a training procedure where, given that the nth learning time X n has duration t, 
the (n + l)th random time is identically distributed to X n conditional on X n < t. In some 
sense, (fl2l) expresses that the information collected at each stage allows the next step of 
the procedure to have a stochastically smaller duration. Alternatively, {X n ,n > 1} may be 
viewed as the sequence of lifetimes of an item that is repaired instantaneously after each 
failure, such that after each repair the duration of the next lifetime is stochastically smaller 
than the previous, due to imperfect repairs and weakening caused by wear. 

From (HU), since T n (x) > 0, for all x > 0 and n > 1, we derive that F n+ i(x) > F n (x) for 
all x > 0. Hence, we conclude that X n > s t X n+ \ for all n > 1. In the following theorem, 
we obtain the same result for a stronger stochastic order. 

Theorem 3 Consider the sequence of random variables {X n ,n > 1} as defined in [T3\) . 
For all n = 1,2,---, we have 


X n >l r X n+ \. 

Proof. From the definition of likelihood ratio order and (1171) . we conclude that 

fn+ iO) 


fn(x) 


= T n (x), n> 1. 


(19) 


The right-hand-side of (|19|) is decreasing in x > 0 and, hence, the claimed result follows. 
Theorem 4 Let n > 1. If X n is ILR then 

X n ^lrj- X n+ 1. 

Proof. From the definition of up-shifted likelihood ratio, it is sufficient to prove that the 
function f n+ i(x + t)/f n (x ) is decreasing in x for all t > 0. First, from (fI71) . we observe that 

fn+ i(s + t) = T + ^f n {x + t) 


fn{x) 

On the other hand, recalling (1151) . we have 


fn(x) 


dx ^ Tn(x + t} 


fn(x + t) 

fn(x) 


< 


[~T n (x + t)f n (x + t) + T n ( X + t)f' n (x + t)\f n (x) - T n (x + t)f n (x + t)f^(x) 

fn(x) 

T n (x + t)[f n {x + t)f n (x) - f n {x + t)f' n {x)\ 
fn(x) 


( 20 ) 


The last Expression in (1201) is negative since, by assumption, X n is ILR, i.e. 

for all 0 < x < t + x. 


fn(x) > fn(x + t) 


fn(x ) fn(x + t)’ 

The proof is thus completed. 

Remark 1 If X is DRHR, i.e. t{x) is decreasing, then T{x) is convex. 
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Theorem 5 Let n > 1. If X n is DLR then X n+ i is DLR. 

Proof. From (fI71) . we need to show that 

log/n+i(®) = log/ n (x) + logT n (x) (21) 

is convex for all x > 0. If X n is DLR then log f n (x) is convex. On the other hand, it is well 
known that if X n is DLR then X n is DRHR and, hence, T n {x) is convex due to Remark |T] 
Since the function log(-) is convex and increasing, logT n (x) is convex. From (121 jl . we have 
that log/„+i(x) is the sum of two convex functions and then the desired result follows. □ 

Remark 2 If X n is ILR then X n+ \ is not necessarily ILR. To show this fact, consider the 
following example. 


Example 4 Suppose that X\ has CDF F(x) = x a , 0 < x < 1 (a > 0). It is easy to see 
that r(x) = a/x, and T{x) = —a logx, 0 < x < 1. It follows that 


f[(x) = a -l 
fi(x) x 


0 < x < 1, 


is decreasing in x G (0,1) for all a > 1 and, hence, X\ is ILR. (Note that T{x) is convex 
according to Remark [lj) On the other hand, we have 


f 2 {x) 1 + (a — l)logx 

f2 (x) xlogX 


0 < x < 1. 


( 22 ) 


It is easy to show that the right-hand side of (|22J) is not decreasing in x for all a > 1 and, 
thus, X 2 is not ILR. 


Proposition 5 Let q(x) be a nonnegative function of x > 0. If q(x)ri(x) is a decreasing 
function of x > 0, then q(x)r n (x) is also a decreasing function of x > 0 for all n = 1,2,_ 


Proof. We just show that under the hypothesis, the function q(x)r 2 (x) is a decreasing 
function of x > 0. From (fT4l) and (fl7l) . we have 


q{x)T 2 (x) = q(x) 


h(x) 

F 2 (x) 


Ti(x)h(x) 

(l + Ti(s))Fi(x) 


q(x)n(x) 


Ti(x) 

1 + Ti{x)' 


(23) 


Since x/{x + 1) is an increasing function of x > 0, and the function T\(x) is decreasing with 
respect to x > 0, then the function Ti(x)/(1 + Ti(x)) is decreasing with respect to x > 0. 
From (1751) . we thus obtain that q{x)r 2 {x) is a decreasing function of x > 0. The rest of the 
proof follows by induction. 


Corollary 1 If X\ is DRHR then X n is DRHR for all n > 2. 

Proof. The proof follows from Proposition [5] by taking q{x) = 1 for all x > 0. 

Let us now consider the following property. 

Definition 3 Let X be an absolutely continuous random variable with support (lx,ux)- 
We say that X has the DRHR in a length-biased sense (LBDRHR) if xt(x) is decreasing 
in x <G (lx,u x )- 
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We remark that a necessary and sufficient condition such that X is LBDRHR has been 
given in terms of stochastic comparison of quantile-based distributions in m- Other results 
on the characterization given in Definition [3] will be the object of a future investigation. 

Corollary 2 If X\ is LBDRHR, then X n is LBDRHR for all n > 2. 

Proof. The proof follows from Proposition [5] by taking q{x) = x for all x > 0. □ 

Consider now the following stochastic order from m- 

Definition 4 Let X and Y be absolutely continuous random variables with reversed hazard 
rates tx(x) and ry(x'), respectively. The random variable X is said to be smaller than Y in 
relative reversed hazard rate order (denoted by X <rrh Y), if ty{x)/tx{x) is an increasing 
function of x. 

For instance, let X and Y denote the lifetimes of two components; given that the 
components have been found to be failed at the same time, then X <rrh Y states that Y 
has been lived longer than X or, equivalently, X aged faster than Y. 

Proposition 6 The sequence of random variables defined in m satisfies X n T rrh X n -y. \ 
for all n > 1. 

Proof. From (fill) . (ITBjl and (1171) . we have 

Tn+l{x) _ fn+l(x)F n (x) _ T n (x) 

T n (x) F n+ i(x)f n (x) 1 + T n (x)' 

We can see that the function T n (x)/( 1 + T n (x)) is decreasing with respect to x > 0, for all 
n > 1, which completes the proof. □ 

We remark that the results stated in Theorem [3] and Proposition [6] are not related to 
each other, since likelihood ratio order does not imply RRH rate order and vice versa. 

Theorem 6 Consider the sequence of random variables {X n ,n > 1} as defined in 113\) with 
the reversed hazard rate functions defined in m, and let Y be an absolutely continuous 
nonnegative random variable with the reversed hazard rate function Ty(x), x > 0. If X i and 
Y satisfy the PRHRM with G{x ) = [Fi(x)] e , x > 0, then X n <rrh Y for all n = 1,2,.... 

Proof. Since t\{x)/ty{x) = 0~ x , x > 0, from © we have that 

= h( x ) = n(x) Ti(x) = x Ti(x) 
t y (x) F 2 (x)ty(x) ty{x) 1 + T\{x) 1 + Ti(x)' 

The function Ti(x)/(1 + T\(x)) is a decreasing with respect to x > 0 and, hence, X 2 <rrh 
Y. The rest of the proof follows by induction. □ 

4 Connection with entropy and covariance 

In this section we obtain some results about the connection between the entropy and the 
inactivity time of the new weighted distribution function considered in the previous section. 
One of the most important measure of uncertainty is the differential entropy introduced by 


11 









Shannon |33j . For an absolutely continuous nonnegative random variable X having PDF 
/, the differential entropy is defined by 

pOO 

H(X) = - f(x) log f(x)dx, (24) 

Jo 

where ‘log’ means natural logarithm and, by convention, OlogO = 0. The entropy H(X) 
gives expected uncertainty contained in f(t) about the predictability of an outcome of the 
random variable X. It is known that in many realistic situations, such as in survival analysis 
and reliability, one has information about the past time, i.e. the time elapsed after failure 
till time t, given that the unit has already failed. The entropy (|24l) applied to conditioned 
random variable is useful to measure uncertainty in such situations. Di Crescenzo and 
Longobardi E3 indeed considered the entropy for the past lifetime, called past entropy at 
time t of X, denoted by 

m = ~L W) log W) ix ' * > 0; 

see also [28]. Furthermore, the concept of dynamic cumulative entropy as an alternative 
measure of uncertainty for the inactivity time was introduced in m and is defined as 

Fix) , Fix), 

( ;> = ~Lm g W) 

1 l' 1 

= -j^jJ^ F {x) l ogF(x)dx-T(t)fl(t), t> 0, (25) 

where fl(-) and T(-) are defined in (JT]) and flU), respectively. Note that 

poo poo 

CS(X) = lim C£[X\ t) = — / F(x) log F(x)dx = / F(x)T(x)dx, (26) 

Jo Jo 

where C£(X) is called cumulative entropy of X. Di Crescenzo and Longobardi [12] also 
showed that the dynamic cumulative entropy and the mean inactivity time are connected 
as 


and, thus, 


C£(X;t) =E\ji(X)\X <t], f > 0 


C£(X)=E[fl(X)]. (27) 

Now, we extend the results of [18] to the case of past time. Based on X\ = X , we consider 
the sequence of random variables {X n ,n > 1}, with the corresponding distributions F n (x) 
as dehned in (HU), and denote by 

Unit) = E[X n | X n <t], t > 0, (28) 


the mean past lifetime of X n , n > 1. Hereafter, we obtain the main results and connections 
between the dynamic cumulative entropy and the reversed hazard rate function. 


Theorem 7 For any t > 0, and for all n = 1,2,..., we have 


T n (X n ) 
7~n (X n ) 


x n <t 


C£(X n ;t) + p, n (t)T n (t), 


cov [X n , T n (X n ) I x n <t\ = T n (t ) [Hn{t) - E(x n )} - C£(X n ;t). (29) 
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Proof. For t > 0, and n > 1, we have 

T n (X, 


IE 


T>{X n 


x n < t 


1 /•* T n (x) 

^n(*) Jo 7Vi(®) 
1 * 


/„(a:)dx 


F ^ j o T n (x)F n (x)dx 
= C£(X n ;t) + p, n (t)T n (t), 

where the last equality is obtained from ([25]) . To prove the second expression, it is easy to 
show that the random variable T n (X n ) is exponentially distributed with unity mean and, 
hence, E[T n (X n )] = 1. Now, consider the following expression for t > 0: 


cov [X n , T n (X n ) I X n <t\ = IE { [X n - E(X n )][T n {X n ) - E(T n (A n ))] \X n <t} 

= E[X n T n {X n ) | X n <t\- n n {t) 

-E(X n )E[T n (X n ) | X n <t}+ E(X n ). 


We have that 


E[T n (X n )\X n <t] = 


Also, we see that 


Fnif) Jo 

= — r 

F n {t) Jo 
Fn+l(t) 
Fnit) 

= 1 + T n (t). 


1 


T n (x)fn(x)dx 


fn+l(x)dx 


E [X n T n (X n ) I X n < t] = —— / xf n+ \{x)dx. 

rn\t) Jo 


Integrating by parts, we can derive 
E[X n T n (X n )\X n <t] = 


1 


F n (t) L 


tF n+ i(t) - [ F n+ i(x)dx 
Jo 


— f [1 + T n (t)\ — 


1 


[1 + T n (x)]F n (x) dx, 


Fnit ) Jo 

and by using (I25|) and after simplification, we obtain 

E[X n T n (X n ) | X n <t\ = fjL n {t)[ 1 + T n (t)] - C£{X n -t). 
Substituting (]3TT) and (]33l) into (1301) . the desired result (1291) finally follows. 

From Theorem [7] the following corollary is derived. 

Corollary 3 Under the conditions of Theorem [?] we have, for n > 1, 

= C£(X n ), 


lim E 

\T n (X n ) 

x n < t 

= E 

\T n (X n )^ 

t—yoo 

_ T~n(X n ) 



_ Tn{X n ) _ 


(30) 


(31) 


(32) 


(33) 

□ 


lirn cov [X n ,T n (X n ) \ X n <t]= cov (X n ,T n (X n )) = -C£{X n ). 

t—yoo 
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Remark 3 Note that the initial random variable was arbitrary selected. Therefore, for any 
absolutely continuous nonnegative random variable X, the following identities hold: 

E (v§j) = C£(X) ’ cov ( X ’ T ( x )) = -C£(X). 

Moreover, for all t > 0, we have 


E 


\T(X) 

X < t 

r{X) 



= C£(X;t) + jl(t)T(t), 


cov[X,T(X)\X <t} = T(t)[p(t) -E(X)] -C£(X;t), 
where fi(t) = E[X | X < t], t > 0, denotes the mean past lifetime of X. 

Remark 4 Consider the sequence of random variables as defined in (11311 . From (1141) . we 
have, for n > 1 , 


F n + iO) = 1 - F n+1 (x) = F n (x) - T n (x)F n (x), x > 0. 

Hence, recalling (1261) . we obtain the following iterative expression for the mean of X n : 

E(X n+1 ) = E(X n ) - C£(X n ), n> 1, (34) 

which also gives a new probabilistic meaning for the cumulative entropy. 

Theorem 8 For all n = 1,2,..., we have 


cov(X n , X n+ i) = cov(X n , X n - jl(X n )) = v&r(X n ) - cov(X n , ft{X n )). ( 35 ) 

Proof. From (P) and (|28|) . we have n(t) = t — and then 

E[X n X n+ i\X n = t] = tE[X n+ i\X n = t] 

= tE[X n \X n <t\ 

= t{t~ fin(t)}, 

for t > 0. Hence, 

E[X n X n+1 ] = E[xZ]-E[X n jl n (X n )\. ( 36 ) 

Moreover from (l34l) and (f36ll . we have 

cov (X n , X n+1 ) = E[X 2 n ] - E[X n jx n {X n )\ - E 2 [X n ] + E[X n ]C£(X n ) 

= var[X n ] - cov(X n ,jl n (X n )) 

— cov(X n , X n )i n (Jf n )). 

The second equality follows from C£(X n ) = E[jl n (X n )\, due to (|27ll . The desired result then 
follows. □ 

Now we use the probabilistic mean value theorem (see M) to obtain an iterative result 
for C£(X n+ 1 ). We first recall the following result (see [12, Equation (12)]). 
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Lemma 1 The derivative of the mean inactivity time of X, given in m, can be expressed 
in terms of the reversed hazard rate function (when existing) as 

pf(t ) = 1 — T(t)p,(t ), t > 0 : F(t) > 0. (37) 


Theorem 9 For the sequence of random variables {X n ,n > 1} defined in \13\) and for all 
n = 1 , 2 ,..., we have 


C£(X n+1 ) = E [fi n+ 1 (X n )} - E \ft' n+ 1 (Z)]C£(X n ), (38) 


where fi' n+1 (t ) can be obtained from [37\ ) and Z is an absolutely continuous nonnegative 
random variable having PDF 


f ( \ _ F n+i{z) ~ F n (z) = F n (z)T n (z ) 
!z{Z ~ E(A n ) - E(X n+ i) “ C£(X n ) 


z > 0. 


(39) 


Proof. Since X n+ \ < st X n and C£{X n+ 1 ) = E[/r n+ i(X n+ i)], the desired result imme¬ 
diately follows from [9, Proposition 3.1 and Theorem 4.1]. Note that the pdf in (1391) is 
obtained from (fTH) and (IM1) . □ 


For the sequence of random variables defined in (1131) one can verify that, for t > 0, 

C£(X n -t) 


Mn+1 (f ) — Tn (f ) T 


i + T n (ty 


n = 1 , 2 ,. 


(40) 


as 


Hence, from (1401) . we can write 

CS(X n+1 ) = C£(X n ) (1 - E \p! n+ 1 {z)]) + E 


C£(X n ; X n ) 

_ 1 + T n (X n ) 


n = 1 , 2 ,. 


5 An integral operator 


Stimulated by some results shown in [18, Section 4], we now define a new operator which 
is dual to the T s operator introduced by Dickson and Hipp [ 8 ]. It is known that the 
Dickson-Hipp operator for any s £ R, denoted by T s f(t), is defined by 

/»oo 

T s f(t) = J e~ s ^f(x) dx, t > 0, 

where f(x) is an integrable function. Applications and properties of T s operator can be 
found in m , 0, and PU, among others. For simplicity, we define an operator T s f(x) for 
an integrable function / and for s S 1, by 


T s /(t) = f e~ s ^f(x)dx, t > 0 . (41) 

Jo 

We can see that, for t > 0, the two operators are related by the following identity: 


T- a f(t) + T 8 f(t) =e st £ s [f], s> 0, 

where 

r oo 

£s[f] = / e~ sx f(x)dx, s > 0 , 

Jo 
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denotes the Laplace transform of the function /. Note that (1411) recalls a generalized Hardy 
operator, similar to that considered in [5, Definition 2], Moreover, T s f(t) can be viewed as 
a convolution-type operator; see [T]. For instance, if s > 0 and f(x ) is a PDF then, due to 
(l4lj) . (l/s)T s /(t) is the convolution between f(x) and an exponential PDF with parameter 

s. 

Suppose that X is an absolutely continuous nonnegative random variable with the CDF 
F. Then the T s operator of F is defined by 

T s F(t) = [ e~ s ^- x) F{x)dx, t > 0, seR. (42) 

Jo 

Integrating by parts, from (14T1) and (H 2 l) . it follows that 

T s f(t) = F(t ) - sT s F(t), t > 0. 


If X is an arbitrary absolutely continuous random variable with the pdf / then, from (1411) . 
we have 

lim e st T J(t) = E[e sX ], 

£—>•00 

Hence, © can be written in terms of the moment generating function of X when t goes 
to infinity. Now, we have an iterative result for T s F n (t). 


Theorem 10 Let {X n ,n > 1} be a sequence of absolutely continuous nonnegative random 
variables with the corresponding CDFs defined in 0 and let T s F n {t ) be defined as in 1 ). 
Then, for t > 0, 


T s F n+ i(t) — [1 + T n (t)] T s F n (t) + T s q n (s, t), n> 1, s € M, (43) 

where q n (s,t ) = T n (t)T s F n (t). 


Proof. From (fl4l) and (pl2p . we have 

T s F n+1 (t) = f e-< t -^F n (x)dx+ [ e~ s ^T n {x)F n (x)dx 
Jo Jo 


T s F n 


T F 

s 1 - n 


(*)+r 

JO 

it) + f 

JO 


0 -s(t-x) 


T n {u)du 


F n (x)dx 


0 -s(t-x) 


/ T n (u)dl 

J X 


F n (x)dx 


+ / e ~ s J~ x) 


T n (u)du 


F n (x)dx 


= [l + T n (t)]T a F n (t)+ f e-'^Tniu) 

Jo 

= [1 + T n (t)]T s F n (t) + [ e- s J- u ^T n (u)T s F n (u)du, 

Jo 

and the proof is completed. 


e -*<M-*)F n (x)dx 


d u 


The following corollary can be obtained from (1271) and ()43|) by setting s = 0. 
Corollary 4 Under the conditions of Theorem 1 1 (A we have, for t > 0, 


T 0 -F n+ i(f) = F n (t)[fi n (t)( 1 + T n (t)) +C£{X n ]t)\, n > 1. 


(44) 
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Figure 1: CDFs of X \, X 2 .... ,X$ when X\ follows the distributions given in Table 1. 

Table 1: Starting distribution functions. 

(a) : F[x) = x/2, 0 < x < 2 

(b) : F(x) = e~ x 2 , x > 0 

(c) : F{x) = e _3 /( ea:_1 ), x > 0 

(d) : F(x) = 1 — e~ x , x>0 

(e) : F(x) = l — e~ x2 , x>0 

(f) : F(x) = 1 — e~ x / 2 (l + x/2 + x 2 /8) , x>0 


5.1 Computational results 

We conclude the paper with few illustrative examples which shed some light on the behavior 
of the sequences of random variables defined in ( 11 .‘ 111 . 

Example 5 Let X\ be uniformly distributed on [0,2]. Then we have T\(t) = — logt/2, 
Ji\(t) = t/2 and C£(Xi;t) = t/ 4, 0 < t < 2 (see |T2]), so that 


E 


( TAX,) 

V r(Xr) 


Xi < t 


C£(X 1 -t) + jl 1 (t)T 1 (t)= t 1 - t - log|, 


and 

covpC, T^Xi) I X 1 <t) = Ti(t)[ni(t) - E(X!)] - C£(Xr,t) = ^ log *- - 
From (|35|) . we obtain 

cov(X 1 ,X 2 ) = cov(X 1 ,X 1 - fii(Xi)) = ^var(Xi) = 


Finally, from (J44J), we obtain 


T 0 F 2 (t) = F 1 {t)[jl 1 (t){l+Ti(t))+C£(X 1 -,t)} = 


°-lpgl 
2 2 


for 0 < t < 2. 
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Table 2: The mean of the starting distribution given in Table 1. 



(a) 

(b) 

(c) 

(d) 

(e) 

(f) 

E(Xd 

1.000 00 

1.772 45 

1.937 91 

1.000 00 

0.886 22 

6.000 00 

E(X 2 ) 

0.500 00 

0.886 62 

1.15166 

0.355 06 

0.506 59 

3.359 35 

E(X 3 ) 

0.180 66 

0.606 80 

0.739 45 

0.104 92 

0.260 35 

1.820 59 

E(X 4 ) 

0.04713 

0.469 75 

0.51103 

0.024 83 

0.118 26 

0.944 89 

E(X 5 ) 

0.009 06 

0.387 74 

0.375 44 

0.004 59 

0.04717 

0.466 17 


Table 3: The cumulative entropy of the starting distribution given in Table 1. 



(a) 

(b) 

(c) 

(d) 

(e) 

(f) 

C£(X{) 

0.500 00 

0.886 23 

0.786 25 

0.644 94 

0.379 63 

2.640 65 

C£{X 2 ) 

0.319 34 

0.279 35 

0.412 21 

0.25014 

0.246 24 

1.538 76 

C£(X 3 ) 

0.133 53 

0.13712 

0.228 42 

0.080 09 

0.142 09 

0.875 70 

C£(X 4 ) 

0.038 07 

0.082 01 

0.135 59 

0.020 24 

0.07109 

0.475 72 


It is difficult to obtain neat analytical results for the sequence of random variables 
{X n ,n > 1} and, therefore, we are forced to proceed via numerical computations. To 
this aim, in Figure 1 we show plots of the cumulative distribution of the random variables 
X\, X 2 ,..., X$ for different starting distribution functions that are listed in Table 1. In the 
figure, the solid line corresponds to the CDF of X\, the large-dashed line corresponds to the 
CDF of X 2 and so on. Moreover, we compute numerically the mean of the recursive random 
variables as well as the corresponding cumulative entropy; see Tables 2 and 3, respectively. 
Recalling (134p . the cumulative entropy is computed as the difference of two consecutive 
means. As expected, the mean of X n decreases when n increases, whereas the cumulative 
entropy decreases when n increases. According to the numerical findings, we expect that 
the cumulative entropies are decreasing for any given starting distribution function. 


Acknowledgements 

We wish to thank and acknowledge partial support from the Ordered and Spatial Data 
Center of Excellence of Ferdowsi University of Mashhad of Iran, from GNCS-INdAM, and 
Regione Campania. Abdolsaeed Toomaj is grateful to the Department of Mathematics, 
Salerno University, for the hospitality during a four month visit in 2014. 

We gratefully thank the anonymous referee for his/her careful reading of the paper and 
for valuable suggestions. 


References 

[1] Abate, J. and Whitt, W. (1996). An operational calculus for probability distributions via Laplace 
transforms. Adv. Appl. Prob. 28, 75-113. 

[2] Bartoszewicz, J. (2009). On a representation of weighted distributions. Stat. Prob. Lett. 79, 1690- 
1694. 

[3] Baxter, L. A. (1982). Reliability applications of the relevation transform. Nav. Res. Logist. Quart. 
29, 323-330. 

[4] Block, H. W., Savits, T. H. and Singh, H. (1998). The reversed hazard rate function. Probab. Eng. 
Inform. Sci. 12, 69-90. 


18 










[5] Bloom, S. (1997). First and second order Opial inequalities. Studia Math. 126, 27-50. 

[6] Burkschat, M. and Navarro, J. (2014). Asymptotic behavior of the hazard rate in systems based 
on sequential order statistics. Metrika, 77, 965-994. 

[7] Cai, J., Feng, R. and Willmot, G. E. (2009). On the expectation of total discounted operating 
costs up to default and its applications. Adv. Appl. Probab. 41, 495-522. 

[8] Dickson, D. C. M. and Hipp, C. (2001). On the time to ruin for Erlang(2) risk process. Insur. Math. 
Econ. 29, 333-344. 

[9] Dl Crescenzo, A. (1999). A probabilistic analogue of the mean value theorem and its applications to 
reliability theory. J. Appl. Probab. 36, 706-719. 

[10] Di Crescenzo, A. (2000). Some results on the proportional reversed hazards model. Stat. Prob. Lett. 
50, 313-321. 

[11] Di Crescenzo, A. and Longobardi, M. (2002). Entropy-based measure of uncertainty in past 
lifetime distributions. J. Appl. Probab. 39, 434-440. 

[12] Di Crescenzo, A. and Longobardi, M. (2009). On cumulative entropies. J. Stat. Plan. Infer. 139, 
4072-4087. 

[13] Di Crescenzo, A., Martinucci, B. and Mulero, J. (2014). A quantile-based probabilistic mean 
value theorem. Probab. Eng. Inform. Sc. 30, 261-280. 

[14] Dimitrov, B., Chukova, S. and Green, D., Jr. (1997). Probability distributions in periodic random 
environment and their applications. SIAM J. Appl. Math. 57, 501-517. 

[15] Gupta, R. C., Gupta, P. L. and Gupta, R. D. (1998). Modeling failure time data by Lehman 
alternatives. Comm. Stat. Theor. Meth. 27, 887-904. 

[16] Gupta, R. C., Kannan, N. and Raychaudhari, A. (1997). Analysis of lognormal survival data. 
Math. Biosci. 139, 103-115. 

[17] Gupta, R. D. and Nanda, A. K. (2001). Some results on (reversed) hazard rate ordering. Comm. 
Stat. Theor. Meth. 30, 2447-2458. 

[18] Kapodistria, S. and Psarrakos, G. (2012). Some extensions of the residual lifetime and its con¬ 
nection to the cumulative residual entropy. Probab. Eng. Inform. Sci. 26, 129-146. 

[19] Kijima, M. and Ohnishi, M. (1999). Stochastic orders and their application in financial optimization. 
Math. Meth. Oper. Res. 50, 351-372. 

[20] Krakowski, M. (1973). The relevation transform and a generalization of the Gamma distribution 
function. Rev. Frangaise Automat. Informat. Recherche Operat. 7, 107-120. 

[21] Lau, K.-S. and Prakasa Rao, B. L. S. (1990). Characterization of the exponential distribution by 
the relevation transform. J. Appl. Prob. 27, 726-729. 

[22] Lau, K.-S. and Prakasa Rao, B. L. S. (1992). Characterization of the exponential distribution by 
the relevation transform. J. Appl. Prob. 29, 1003-1004. 

[23] Lehman, E. L. (1953). The power of rank tests. Ann. Math. Statist. 24, 28-43. 

[24] Li, S. and Garrido, J. (2004). On ruin for the Erlang(n) risk process. Insur. Math. Econ. 34, 391-408. 

[25] Li, Y., Yu, L. and Hu, T. (2012). Probability inequalities for weighted distributions. J. Statist. Plann. 
Infer. 142, 1272-1278. 

[26] Mudholkar, G. S. and Hutson, A. D. (1996). The exponentiated Weibull family: some properties 
and a flood data application. Comm. Stat. Theor. Meth. 25, 3059-3083. 

[27] Mudholkar, G. S., Srivastava, D. K. and Freimer, M. (1995). The exponentiated Weibull family: 
a reanalysis of the bus-motor-failure data. Technometrics 37, 436-445. 


19 


[28] Muliere, P., Parmigiani, G. and Polson, N. (1993). A note on the residual entropy function. 
Probab. Eng. Inform. Sci. 7, 413-420. 

[29] Navarro, J., Sunoj, S. M. and Linu, M. N. (2011). Characterizations of bivariate models using 
dynamic Kullback-Leibler discrimination measures. Stat. Prob. Lett. 81, 1594-1598. 

[30] Psarrakos, G. AND Navarro, J. (2013). Generalized cumulative residual entropy and record values. 
Metrika 27, 623-640. 

[31] Rezaei, M., Gholizadeh, B. and Izadkhah, S. (2015) On relative reversed hazard rate order. 
Comm. Stat. Theor. Meth. 44, 300-308. 

[32] Shared, M. and Shanthikumar, J. G. (2007). Stochastic Orders. Springer, New York. 

[33] Shannon, C. E. (1948). A mathematical theory of communication. Bell. Syst. Tech. J. 27, 379-423, 
623-656. 


20 


